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Abstract: Motivated by the gluon condensate in QCD we study the perturbative ex- 
pansion of a gauge theory in the presence of gauge bosons of vanishing momentum, in 
Qh. the specific case of an abelian theory. The background is characterised by a dimensionful 

parameter A affecting only the on-shell prescription of the free (abelian) gluon propagator. 
When summed to all orders in gA the modification is equivalent to evaluating standard 
Green functions in a pure gauge field with an imaginary gauge parameter oc A. We show 
how to calculate the corresponding dressed Green functions, which are Poincare and gauge 



covariant. We evaluate the expressions for the dressed quark and qq propagators, im- 
posing as boundary condition that they approach the standard perturbative form in the 
short-distance limit {\p 2 \ — ► oo). The on-shell (p 2 = m 2 ) pole of the free quark propagator 
is removed for any A > 0, and replaced by a discontinuity which vanishes exponentially 
with p 2 . The dressing introduces an effective interaction between quarks and antiquarks 
which is enhanced at low relative 3-momentum. Further study should allow to identify the 
(bound) eigenstates of propagation and determine whether they define a unitary 5-matrix. 
When the quark mass is zero there is a euclidean propagator solution which breaks chiral 
symmetry spontaneously. We study some aspects of the massless pion contribution to axial 
vector correlators and derive the 7rqq form factor. 
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1. Introduction 



The perturbative expansions of field theories such as QED and QCD are defined by the 
lagrangian once the boundary conditions are specified. It is customary to use the pertur- 
bative vacuum |0) as the state around which the expansion is made. This state is empty, 
as by definition dfe|0) = for any annihilation operator a^. \in) and \out) states with 
particles are built by applying creation operators on the perturbative vacuum: aj[,|0), etc. 

The perturbative expansion works quantitatively for QED to an amazing precision, 
indicating that the true QED ground state is quite close to the perturbative vacuum. In 
addition to providing quantitatively reliable results for QED, the perturbative expansion 
more generally provides analytic and unitary scattering amplitudes - as expected for S- 
matrix elements in a physical theory. 

The perturbative expansion of QCD is qualitatively similar to that of QED. However, 
we know from experiment that quarks and gluons (unlike electrons and photons) do not 
appear as asymptotic scattering states. The S-matrix of QCD formulated in terms of 
hadron states should be analytic and unitary. QCD Green functions of quarks and gluons, 
on the other hand, are not simply related to the S-matrix and their analyticity and unitarity 
properties are largely unknown. 

The failure of perturbative QCD to describe the observed hadrons may be caused by 
the QCD ground state being a quark and gluon condensate. An expansion around the 
perturbative vacuum would then be inadequate to describe the long distance propagation 
of quarks and gluons. This need not imply that the coupling a s is large - in fact, there are 
plausible arguments Q that the coupling freezes at a moderate value in the long distance 
regime. The situation may be compared to the propagation of a high energy electron in 
matter: multiple scattering limits its range even though a em is small. 

It is thus interesting to contemplate perturbative expansions around states other than 
the perturbative vacuum, which might be closer to the ground state of QCD. Such expan- 
sions could indicate how an analytic and unitary S"-matrix is obtained when the asymptotic 
states are not simply given by the fields of the lagrangian. It is of course desirable that the 
boundary condition does not break basic symmetries such as Poincare and gauge invariance. 

The Feynman ie prescription for free quark and gluon propagators is a consequence of 
expanding around the perturbative vacuum. This may be illustrated by the expression of 
the free scalar Feynman propagator when the ground state contains a particle of momentum 



k, 






which in momentum space is 




(1.2) 



+ 



(27T) 4 



[5 3 (p-k)5(p° -E k ) + 5 3 



(p + k)5(p° + E k )] 



N 
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The term in the second line arises from interference between the background and propa- 
gating scalars for p = ±k. Only the on-shell propagator is affected, since the background 
particle is on-shell by definition. 

In this paper we study the effect on abelian gauge theory of using an (abelian) gluon 
propagator of the form 



D^{p) = -g 



fJbV 



jr + le 



Cl-3) 



The modification oc A 2 is analogous to the one in ( |1.2| ) for k = 0. Thus the last term 
in ( |1.3|) signals additional zero-momentum gluons in the \in) and \out) states. We shall 
not here be concerned with the precise boundary condition implied by this modification 
of the gluon propagator (see the appendix of for a discussion of this). It turns out, 
however, that the effect of the modification is formally equivalent to that of a pure gauge 
field background with an imaginary gauge parameter (see Eq. (|2.9| ) below). 

In ||, Q we showed that the modification ([0^) can be obtained by assuming a constant 
background field and averaging over all components of this field with a gaussian weight. It 
was possible to calculate exactly the dressed quark and gluon propagators in the limit of 
a large number of colors (N c — > oo), i.e., for all planar Feynman diagrams. The present 
study is restricted to the abelian case, and we use a different method to sum the large 
number of (planar and non-planar) diagrams. 

Modifying the (on-shell) gluon propagator only at vanishing momentum (p = 0) pre- 
serves Poincare invariance and enables us to evaluate exactly the effects of the modification. 
Thus we get a novel perturbative expansion, with Green functions at each order in a dressed 
to all orders in gA. We impose as boundary condition that the Green functions reduce to 
their standard perturbative expressions at short distance (p 2 — ► ±oo). At long distance 
\p 2 \ < (gA.) 2 they turn out to have a qualitatively different (exponential) behavior. 

We view this study as a possible step towards understanding the long-distance proper- 
ties of QCD. Even though our abelian example is at best a toy model, the explicit results 
we obtain for the perturbative expansion with a non-trivial boundary condition may be 
useful for understanding properties of gauge theories in a non-perturbative sector. The 
vacuum of a perturbative expansion should be such that the essential physics is manifest 
already at low orders of the expansion. We do not know whether the present approach can 
fulfill this criterion when applied to QCD. Phenomenological studies || using modifications 
of the gluon propagator similar to Ql.3| ) are, however, encouraging in this regard. 

In section 5 we note that for vanishing quark mass there is a quark propagator solution 
which breaks chiral symmetry (XSB). The symmetry breaking term vanishes exponentially 
for p 2 — > — oo. Hence the XSB propagator gives well-defined results in Wick-rotated (eu- 
clidean) loop integrals. Using the methods developed in this paper we verify the axial 
Ward identity, study the massless pion contribution to axial vector correlators and derive 
the 7rqq form factor. 

2. Dressed quark propagator 

In this section we derive the (abelian) quark propagator in the dressed tree approximation 
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defined in Fig. 1. We sum all (abelian) gluon loops but include only the last term oc A 2 
in each gluon propagator (|1.3|), As this part of the propagator carries no momentum the 
loop integrals are eliminated and the sum can be exactly evaluated. 



+ 



+ 



Figure 1: The dressed quark propagator to all orders in gA (Ihs.) is given by the sum of all gluon 
loop corrections (dashed lines). Only the modification oc A 2 is included in each gluon propagator 



(1.3). Self-interactions (via quark loops) of zero- momentum lines are neglected. 



Using standard Feynman rules, the sum shown in Fig. 1 for the dressed tree propagator 
begins as (an ie prescription is implied) 

i 2i{gK) 2 ( , 2m 2 



S(p) 



■p 1 — m (p 2 



m 



2\2 



p 1 — m 



(2.1) 



Since the parameter A will always occur in the combination gA it is convenient to define 



fx = gA. 



(2.2) 



Furthermore, for dimensional reasons the 
propagator (and any Green function) can 
be written as a function of the dimension- 
less ratios p//i and m//x, times the power 
of /j, corresponding to its dimension. We 
therefore change our notation as follows: 

p m 

>P , — 

\i n 

In the following, momenta and masses will 
be always dimensionless, unless otherwise 





a p 



Si (p) 



Si-i(p) 



m 



(2.3) 



Figure 2: Recurrence relation (2.6) for the quark 
propagator dressed with I gluon loops. The gluon 
line starting at the first vertex on the right is to 
be inserted at all positions between the remaining 
vertices. 



indicated. With this new notation the dressed tree quark propagator can be expressed as 

1 i 



m 



(2.4) 



where Se(p) denotes the £doop contribution to S(p). 

For I > 1 Se(p) satisfies the recurrence relation shown in Fig. 2, where the arrow on 
the (dashed) gluon propagator indicates that it should be inserted at all possible locations 
in Si-i(p). At each location this turns a free propagator Sq(p) into Sq(p)(— ig^f a )So(p)- 
According to the identity 

d 



g_ 

fM dp a 



S {p) = S (p) {-igja) S Q (p) 



(2.5) 
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the same is achieved by differentiating So(p). By differentiating S^_i(p) all insertions of 
the gluon propagator are accounted for. Explicitly, 



St(p) = ~ T 1 — (-^7 a )("A 2 ) {- 9 -\ S^{p) 



(2.6) 



fx p 1 — m v " 1 /s / \ fi J dp a 4 lvr/ p 1 — rn 

where fi = ^y a d/dp a . From ( |2.4j ) and ( |2.6| ) we find that S(p) satisfies the differential 
equation 



(j/> — m 



S(p) 



(2.7) 



This is actually the Dyson- Schwinger equation of Fig. 3, where the dressed qqg vertex is 
obtained by differentiating the dressed quark propagator. 




+ 




Figure 3: Dyson-Schwinger equation for the dressed quark propagator. 

It is interesting to compare the condition ([2]?]) (in coordinate space) with the equation 
of motion for the propagator in a background field A v , which reads 1 



m — —4 (%) 



S{x) = ifi 3 5 4 (x) 



The two previous equations are equivalent for 



gA v {x) = i[ix v = ind v (\x l ) 



(2. 



(2- 



which formally is a pure gauge with an imaginary gauge parameter. The factor i follows 
from our modification of the on-shell gluon propagator in ( |1.3| ), and implies that the field 
Q2.9| ) actually influences the physics. The "gauge transformed" Green functions differ from 
the free ones not just by a phase but in absolute size 2 . As a consequence, a general solution 
of the differential equation (|2.7|) (which we shall also call the equation of motion) for the 
dressed quark propagator grows exponentially at large p 2 . We select a physical solution by 
imposing as boundary condition the perturbative result at 1/p 2 — > 0, i.e., at large p 2 in all 
directions on the physical sheet. 



The equation of motion (2.7) is transformed into two coupled ordinary differential 
equations by using Lorentz invariance to express S(p) as 



S(p) = a{p 2 )i> + b(p 2 ) . 
For Rep 2 > the physical solution is given by the integral representation 

,2 



S( P ) 



2/i 



(j/> + m 



dt exp 



o 



t 



m 

T+t 



(2.10) 



(2.11) 



1 Recall that we are using dimensionless coordinates, fix — > x. 

2 This is reminiscent of the original definition by Hermann Weyl in f 9f 8 of gauge symmetry as invariance 
under a change of scale. 
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We may check that S(p) satisfies the equation of motion (| 
O p = (j/> - m - p + m - p ) = (p 



is a scalar in Dirac space and that 



Op exp 



2^ 



m 

T+t 



d 



by noting that the operator 
m 2 (2.12) 



2-<j(l + t) 2 exp 



2^ 



m 

T+t 



(2.13) 



When inserted into the equation of motion ( |2.7| ) the integrand of the propagator (Ernp thus 
turns into an exact derivative. The integral is then trivially evaluated and (|2.7D verified. 



It is readily seen from fl2.ll ) that for Rep 2 > m 2 we recover the standard free quark 



propagator, either in the p - 
limit at fixed p 2 and m 2 : 

S(p) 



limit at fixed (dimensionful) p and m , or in the p — > oo 



S(p) 



% 

i 



(2.14) 



^i^o f — m p' 

The dressed propagator is regular at p 2 = m 2 for any p > 0. The residue of the pole of the 
free propagator thus vanishes discontinuously with the introduction of a background. 

Propagator representations which converge in the upper (Imp 2 > 0) and lower (Imp 2 < 
0) half planes, respectively, may be found by rotating the integration contour in (|2.11| ) by 
T90°, 



S(p) 



±2> 



m 



dt exp 



m 



V 



IT it 



(2.15) 



From these expressions it is apparent that S(p) approaches the free propagator for |p 2 
in any direction on the p 2 plane. 

If the integration contour of the representation ( 2.1 lj ) is rotated by —180° the t- 
integration will range from to — oo and converge for Rep 2 < 0. The path of integration 
passes below the essential singularity of the integrand at t = — 1. This expression is an ana- 
lytic continuation of S(p) where p 2 passes from positive to negative values above the origin. 
Similarly, if the contour of ( p. 11 ) is rotated by +180° the contour passes above t = — 1. 
The discontinuity of S(p) for p 2 < is given by the difference of the two expressions, 

I r r + / ™2 

Disc S(p 2 < 0) = -— (jt> + m-p) I (it exp 
2p Jc 

where C is a closed contour circling the essential singularity at t 
(counter-clockwise) direction. 

The integral is evaluated by extracting the simple pole contribution at t = — 1, whose 
residue turns out to be a Bessel function, 



m 

i+1 



(2.16) 

1 in the positive 



Disc S(p 2 < 0) 



7T 



/i 



(P 1 



m 




Ji (mV-p 2 ) exp [i(p 2 




(2.17) 



The discontinuity vanishes exponentially for p 2 — > — oo, in agreement with the perturbative 
propagator having no discontinuity. Regarded as an analytic function of p 2 the disconti- 
nuity is seen to increase exponentially for p 2 — ► +oo. This shows that the propagator is 
exponentially behaved on sheets other than the physical one. 
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For p 2 —* the dressed propagator ( [2.11 ) behaves as 



S(p 2 ) 



ie™ 2 ' 2 



^ + ^o gP 2 + 0{{ P 2 f) 



(2.18) 



Thus the leading singularity oc ft/p 4 and is exponentially enhanced in m 2 . We recall that 
p and m are dimensionless - thus the above limit corresponds to p 2 <C [J? in terms of 
dimensionful momentum. 



For a massless quark the propagator ( |2.11[) reduces to 



S{p) 



m=0 



p up 2 



% ( 2 

—7 (1+^7 



(2.19) 



Thus there is only a first order (in /j?) correction to the massless free propagator. We may 
verify that for m = the recurrence relation (2.6) indeed terminates for p ^ 0: 



S*(P) = Tin. 



1 

J 



2%i> 
hp 4 



l —T 5 (P) 

Hp 



(2.20) 



To conclude this section, it is interesting to note that the scalar integral in the repre- 
sentation ( p. 11 ) of the quark propagator actually corresponds to the quark propagator as 
evaluated in scalar abelian gauge theory. We show this in Appendix A, see (A.4). If one 
recalls the background field interpretation (see ( ]2i 
cause for a vanishing field strength tensor (ift — ^.A 



and (2.9)) this is understandable be- 



m)(ift — +m) 



-m 



so that the scalar equation of motion follows from the Dirac equation. The expansion of the 
scalar integral in ( |2.11 ) in powers of /i 2 is given in ( |A.6| ) and is asymptotic - the expansion 
coefficients grow factorially and the series thus diverges for all values of \i. 



3. Alternative derivation of the dressed quark propagator 

In the previous section we demonstrated that the dressed quark propagator satisfies the 



differential equation (|2.7|). Making use of its Lorentz structure ( 2.10 ) we could solve the 



equation and find the expression ( [2. 11 ) for the propagator. As we shall see in the next 



section, higher point Green functions satisfy similar equations of motion, but depend on 
several invariants and have a more involved Lorentz structure. Thus it seems difficult to 
solve the differential equation even in the case of a double fermion propagator. 

Here we present an alternative derivation of the dressed quark propagator, which is 
easier to apply in the case of higher point Green functions. We make use of the observation 
that the differential equation ( |2.7D is equivalent to the equation of motion of the propagator 
in a background field ( |2.8| ) which is a pure gauge ( |2.9j ), albeit with a complex gauge 
parameter. This allows us to immediately write down a solution of the differential equation. 
While this solution grows exponentially with p 2 , we find the physical solution by adding a 
solution of the homogeneous equation 

y-m-p p )S H (p) = ■ (3.1) 
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The equation of motion denned by (2J3) and (|2.9|) in coordinate space, 



(ip x -m- ifi)S{x) = in 3 5 4 (x) , (3.2) 

is satisfied by 

S(x) = exp(ix 2 )S'o(x) , (3.3) 

where the free propagator obeys 

(ip x -m)S (x) =in 3 5 4 (x) . (3.4) 

The propagator S(x) grows exponentially with x 2 and thus cannot be Fourier transformed 
to momentum space. In order to have well-defined expressions we temporarily use euclidean 
metric and switch back to Minkowski space at the end. Our aim is to find a well-behaved 
solution of the equation of motion in Minkowski space - thus the present excursion to 
euclidean metric may be regarded as a purely technical device. 

Making the zeroth component of minkowskian 4-vectors and 7-matrices imaginary 
relates their dot products to the euclidean equivalents through substitutions of the form 

p 2 - -p 2 , i> - -i> , P p - P p , fr". Y} = , 75 - 75 • (3.5) 

We shall flag all equations below where the euclidean metric is used by 'eucl. '. Thus the 
equation of motion 

{ip x — m + ifc)S(x) = ifi 3 5 A (x) eucl. (3-6) 

is solved by the momentum space propagator 

* f A i„-x i~ m 



Snip) = J d 4 xe-^-exp(-Ix 2 )^ J 4^^" 



V J 



fi J (27r) 4 q 2 + m 2 

4r\2 ex p[~i(g-p) 2 ] 2 ™2 eucL ( 3 - 7 ) 

This dressed euclidean propagator is a gaussian averaged version of the free one 3 - with the 
smearing scale set by our parameter (i which defines the dimensionless momenta according 
to d2.3| ). Already from this expression we can see that the propagator will have a bad 
asymptotic behavior at large p 2 in Minkowski metric. Namely, ( |3.8| ) defines S(p) as a 
function of p, which we are free to analytically continue to negative p 2 . While the integral 
always converges, the factor exp(— \p 2 ) in the integrand implies that the propagator grows 
exponentially for p 2 — > —00 {eucl.). The subscript 'D' on the propagator fl3.8|) indicates 
that this propagator is distinct from our previous solution ( |2.11 ). 



^In the notation of our previous work would be the background vector potential 
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Introducing the Feynman parametrisation 



1 



q 1 + m 2 



du exp 



u 



{q 2 + m 2 ) 



eucl. 



and doing the integral over q gives 



i f°° du 



u 

2 V 1 + u 



+ m 



eucl. 



(3. 



(3.10) 



Finally, we change the integration variable to t = —u/(l + u) and return to Minkowski 



metric using the translation (3.5). This gives 

s d(p) = 7r(p , + m -h) / 

Jo 



dt exp 



'2 [P 



m 



1 + t 



(3.11) 



The only difference with our solution ( 2.1 1] ) is in the upper integration limit. The identity 



(|2.13[) ensures that also Sd(p) satisfies the equation of motion (|2.7|) , since the integrand 
vanishes at the endpoint t = —1, just as it does at t = oo for S(p). However, Sr>(p) grows 
exponentially for p 2 — > oo. 



The vanishing of the integrand of ( 3.11 ) at t = — 1 and t = oo (for p > 0) implies that 



Sh(p) = ^(^ + m 



dt exp 



m 



l + t 



(3.12) 



satisfies the homogenous equation (|3.1|) . Our physical solution (|2.11| ) is given by S(p) = 
Sd(p) + Snip)- As we demonstrated in the previous section, S(p) approaches the free 
propagator as \p 2 \ — ► oo in all directions on the physical sheet. 



4. Dressed qq propagator 

We next consider the dressed double propagator of a qq pair (Fig. 4) . Without dressing it 
is just a direct product of two free propagators, 



iGf > p °(k,k) 



i ft + m 



/i 



k 2 



77T 



a/3 



i ft + m 
fi k 2 — m? 



pa 



(4.1) 



where the subscripts on the double propagator G n ,n(^j k) indicate that the number of gluon 
vertices on the q line (of momentum k) is n and on the q line (of momentum —k) is n. In 
order to avoid spelling out the Dirac indices we shall write the direct product in brackets, 
thus 

iG ,o = [S (k)} [S (k)} . (4.2) 
We wish to evaluate the fully dressed qq propagator 



(4.3) 



n=0n=0 



imposing as boundary condition that it approach the free result (|4. l[) in the short distance 
limit. Due to gluon exchanges between the q and q lines G(k, k) will not factorize into 



-8- 



a product of propagators as is the case for Go,o in (| 
exchanges we have, e.g., 

iG 0lii = [S (k)] [S n (k)} 



Only in the absence of such 

(4.4) 



where Sn(k) is the contribution to S(k) which has n vertices on the q line. 

G n! n{k,k) satisfies a recurrence relation which is analogous to that of the single pro- 
pagator in Fig. 2. In diagrams with at least one vertex on the q line (n > 1) the gluon 
propagator from the rightmost vertex must end either to the left on the same line, or on 
the q line (Fig. 4). According to ( |2,5| ) all such insertions are obtained by differentiating 
the double propagator with two vertices less 4 : 



iG n n (k, k) 



fx Jt 



m 



-H<n^)(-A 2 ) 



d 

~dky 



d 

iG n - 2 ,n{k, k) + ^jG n ^i iS _i(fc, k) 



(4.5) 



m 



\G n - 2fi (k, k) + ft- k G n - h ,- n -i{k, k)) (n > 1) 



We may now sum both sides of (|4.5|) over n and n, making use of (^4.4[), 



12 G nAk, k) = G(k, k) + i[S (k)] [S(k)} 

n=l n=0 

Here we assumed the momenta k and k to be 
independent. We may instead regard k as a depen- 
dent variable by keeping the total momentum p of 
the qq system, 

p = k-k (4.7) 

fixed when differentiating wrt. k. Then <j) k will act 
also on A: = k—p and we get the equation of motion 
(or Dyson-Schwinger equation) for G in the form, 

(jt-m- ft k )aa'G a >/3, P a(k, k) = — 8 a p S pa (k) (4.8) 



1 



m 



\ + i/)- k )G{k,k) . (4.6) 




k G(k,k) k 

Figure 4: Rhs. of the recurrence 



relation (4.5) for the qq propagator. 
Analogous to the rhs. of Fig. 2. 



where we exhibit the Dirac structure for clarity. 

5 fc is (in an operator sense) the inverse of the dressed 



According to (2/7), jt — m 
propagator S(k). In (|4.Sj ) it similarly eliminates the quark line from the qq propagator 
G. This is non-trivial since (as we shall see below) the dependence of G(k, k) on k and k 
cannot be factorized. 

Equations analogous to ( |4.8j ) where the inverse dressed propagator acts on the other 
Dirac indices of G can be similarly derived. However, since G(k, k) depends on three 
invariants (which may be taken to be p 2 ,k 2 and k 2 ), and many Lorentz structures with 
four Dirac indices can be built from the two independent momenta p and k, it appears 
difficult to solve the equations of motion directly. We shall instead turn to the alternative 
method demonstrated in the previous section. 



4 Terms with negative indices are understood to vanish. 
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Let us first check that the double propagator Go(k, k) analogous to (3.7) is a solution 
of the equation of motion (JD1) in euclidean metric (cf. ( |3.5[ )). 



G D (k,k) = — i ( . 

u 2 J (2vr) 



d 4 Q , i 2n 
2-exp(-ig ) 



+ a + m 



+ + m 



eucl. 



(4.9) 



As in ( |4.2j ) brackets are used to indicate Dirac structures. The derivative appearing in 
flOj) gives 



p k G D {k,k) 



d 4 q 
u 2 J (2vr) 2 

d 4 q 



1 



(27T) 



2 exp(-ig 2 ) 



7v- 



1 



+ + m 



1 



^ + + m 



(4.10) 



+ + m 



eucl. 



In the first line the derivative d k (which operates equally on k and k) can be replaced by 
d q . The result on the second line then follows after a partial integration. Thus (4.8) is 
indeed satisfied, 



+ m + p k )G D (k,k) = -S D (k) 



eucl. 



(4.11) 



with the antiquark propagator on the rhs. being the euclidean one (3.7). We anticipate 
that Go will not be a solution with good short-distance behavior in Minkowski space, but 
must be corrected by a solution of the homogeneous equation of motion. 

It is helpful to extract the Dirac structure from the integrand as in ( |3.8[) . Similarly to 
(4.10) we note that 



d A q 
a 2 J (2vr) 2 

d 4 q 



(k + q) 2 + m 2 



1 



u 2 J (2vr) 2 



exp(-±<? 2 ^ 



(k + q) 2 + m 2 



1 



(4.12) 



f, + i + m 



eucl. 



It follows that 
G D {k,k) = 

~ l'~ 



d 4 q eM-W) 
(27r) 2 (k + q) 2 + m 2 

d 4 q 



+ i + m 



exp(-i<? 2 ) 



(2tt) 2 {(k + q) 2 + m 2 }{(k + q) 2 + m 2 } 
m + p k ]G DI (k 2 ,k 2 ,p 2 ) eucl. (4.13) 



From the derivation it is clear that ft k in the first bracket acts on the $ in the second 
bracket as well as on Gdi- 
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Introducing Feynman parameters as in ( p.9|) and doing the ^-integrals gives 



G DI {k\k z y) = - 



du dv 



(1+u + v) 



uk 2 + vk 2 + uvp 2 1 



exp ; ^ (u + v)m" eucl. 

2(l + u + v) 2 V ' 

(4.14) 

This expression may now be converted back to Minkowski metric using (|3.5| ). At the same 
time we adjust the integration range as we did for the quark propagator ( |3.11 ). We note 
that the integrand in (4.14) vanishes both for u, v = oo and for u + v = —1. Changing the 
integration range as 



/>CO f*— 1 

/ dudv — ► / dudv 6(1 + u + v) 
Jo Jo 



(4.15) 



will therefore amount to adding a solution of the homogenous equation of motion. Changing 
also the signs of u and v we arrive at the Minkowski expression 



G(k, k) = — [^ + m-p k ][^ + m- p k \ G^k 2 , ~k 2 ,p 2 ' 
f 1 



(4.16) 



with 



Gi(k 2 ,k 2 ,p 2 ) = - dudv 



6(1- u- v) 
(1 — u — V ) 2 



exp 



uk 2 + vk 2 — uvp 2 



1 



1 — (u + v)m 2 

2(1 -u-v) 2 K ' 



(4.17) 

The expressions ( [4.16 ) and (4.17) may be regarded as Ansdtze for the dressed qq 
propagator, motivated by the analysis in euclidean space. We next verify that the equation 
of motion Q4.8| ) is indeed satisfied. Multiplying G(k, k) by ft — m — p k we find the operator 
Cfc defined in ( |2.12j ). Since commutes with + m — ft k , 

[O k j + m-p k ]=0 (4.18) 

it acts only on Gi(k 2 ,k 2 ,p 2 ). The identity 



(O k + 2d u ) 



9(1- u-v) 
(1 — u — v) 2 



exp 



uk 2 + vk 2 



uvp 2 1 , 2 
+ + v ' m 







(4.19) 



2(1 -u-v) 

ensures that O k turns the integrand of Gi into a derivative of u, so that the li-integral is 
evaluated by substitution. The integral over v assumes the form of the single propagator 
( 2.11[ ) through the change of variable t = v/(l — v), showing that fl4.16j ) is a solution of 
Q4.8D - We note that ( |4.17| ) is well-defined provided uk 2 + vk 2 — uvp 2 > when 1 — u — v — > 0. 
This requires k 2 , k 2 > and p 2 < (\fW + a/P) 2 . 

The equation of motion ( |4.8| ) ensures the Ward-Takahashi identity for the (non-am- 
putated) gqq vertex function formed by contracting the /3 and p indices in Fig. 4 with 



r 



Ka(k,k) = 



+ m- P k ) a /3l0p(^ + m- p k ) p(J Gi(k 2 , k 2 ,p 2 ^ 



(4.20) 



Thus, using (| 



p u T u (k,k) = -= (ft + m 
H 2 



-S(k) - -S(k) 
u n 



m 



W-^-m-p^^ + m- p k ) Gj(k 2 , k 2 ,p 2 ) 

(4.21) 



- 11 - 



follows using (|4.8| ) and the analogous equation of motion with the — m — ft k operator. 
More generally, we expect our dressing to preserve gauge symmetry since it corresponds to 
adding on-shell external (abelian) gluons, whose helicities are contracted using the metric 
tensor. 

Similarly, we may verify the axial Ward identity 



for the axial vector vertex 



-2imT 5 (k, k) + -S(k) l5 + -j 5 S(k) 



where 



fi 2 



T 5 (k,k) = -4r(^ + m 
A* 



(4.22) 



(4.23) 



(4.24) 



Taking the short distance limit k 2 — ► oo and k 2 —* oo in Q4.17| ) we have u < 1/k 2 , 
1/k 2 and the integral can be evaluated to give 



iG(k,k) 









'% 


— *oo 

fe^ — ►oo 


_/x k 2 




H k 2 



(4.25) 



To prove the corresponding relation for k 2 ,k 2 — ► — oo would require extending the represen- 
tation analytically. Here we shall only consider the expression for G(k, k) in the massless 
case, m = 0. Replacing the integration variable v in (4.17) by w = 1/(1 — u — v) the integral 
over w can be done exactly. A further change of variable to t = 2uk 2 /(l — u) then gives 



G!(k 2 ,k 2 ,p 2 



m=0 



exp(-t/4) 



t 2 + ttk ■ k + ik 2 k 2 
exp(-t/4) 



o (t + 2k-k) 2 - \(k 2 ,k 2 ,p 2 ) 



(4.26) 
(4.27) 



where A(a, 6, c) = a 2 + 6 2 + c 2 — 2a6 — 2ac — 26c. The integral converges for all values of the 
invariants (given a prescription at the singularities of the integrand) and the short distance 
limit ( 4.25| ) is seen to follow from ( |4.26j ) for k 2 , k 2 growing in any direction of their complex 
planes. We recall that X(k 2 , k 2 ,p 2 ) is related to the relative 3-momentum of the qq pair, 

X(k 2 , ~k 2 ,p 2 ) =A[(k- p) 2 - k 2 p 2 ] = 4(p ) 2 k 2 , (4.28) 

where the second equality holds in the rest frame of the qq pair, with p = (p°,0). 



5. Spontaneous chiral symmetry breaking 

The propagator solutions which we discussed so far maintain chiral symmetry for vanishing 
quark mass, m = 0. Here we consider m = solutions which spontaneously break the chiral 
symmetry of the lagrangian (XSB). Because of the presence of dimensionful constants (such 
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as f w and ip n , see ( 5.18| ) and ( 5 . 1S| ) ) in the following discussion, for clarity we consider in 

b( P 2 ) 



this section only dimensionful momenta (i.e., not scaled by [i as in (2.3)). 

The differential equation (|2.7D for the quark propagator S(p) = a{p 2 )f 
m = 



y-IM 2 p v )(af + b) 



implying 



p 2 {a 



2/iV) 



i>{b 



2u 2 b>) 



4/x 2 a 




where a' = da(p 2 )/dp 2 , and similarly for b(p 2 ). 

The solution which approaches the standard perturbative propagator for \p 2 \ 
given by (12.191 ), 

2^ 

i2 



a(p 2 



1 + 



P 



b(p 2 







The solution of the homogeneous equation fl5,2|) for b(p 



S x (p) = b(p 2 



Aexp 



2 J X 



is for 



(5.1) 



(5.2) 



OO IS 



(5.3) 



(5.4) 



where f x ^ is an arbitrary constant (with the dimension of mass), breaks chiral symmetry 
and is exponentially suppressed at short distance, p 2 — > — oo. It is not well behaved for 
Rep 2 — ► +oo and was therefore excluded in our previous analysis. However, S x (p) is 
meaningful in euclidean space, i.e., in Wick-rotated loop integrals. Keeping this in mind 
we next consider some implications of the dressed propagator solution (with m = 0), 

.2» 



where S x (p) is given by 



S(p) 
~M and 



( P '- fI %)S I (p) + S x (p) 



Slip) = -o 



The value of the quark condensate is 

d 4 p 



<0|qq|0> 



(2iry 



Tr Sip) = 4 



fx 
/i 2 



d 4 p 
(2^Y 



cxp 



P 

2ii 2 



7T Z 



(5.5) 



(5.6) 



(5.7) 



ip . 



where in the second equality the integral is Wick rotated to euclidean space, p° 

It is crucial for an approximation scheme to maintain the axial ward identity (AWI) 

( 4.22 ) H in order to correctly describe the physics of XSB JjJ. For gauge theories one such 

expansion is previously known, based on a systematic truncation of the Dyson-Schwinger 

equations || . Thus we first verify that the AWI is valid in our dressed perturbation theory 

(at the Born order in which we are working). 

We already noted above that the axial vertex ( 4.23 ), corresponding to the propagator 

solution with f x = 0, does satisfy the AWI. In order to construct the axial vertex with 

f x ^ (but m = 0) we note that in analogy to ( |3.7| ) 



Sip) 



1 

7 l 



f d*q 


(q-p) 2 ' 


1 (2n) 2 6XP 


2/i 2 



So(q) 



eucl. 



(5.8) 
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for the choice of undressed propagator 



So(q) 



1 



- + ^/ x (27r)^ 4 (g) 



eucl. 



(5.9) 



where eucl. denotes the euclidean metric ( |3.5| ). The dressed qq propagator is now obtained 
by adding the XSB contribution to the bare propagators in the expression ( fO| ) of Gn(k, k), 

G D (k,k) + G Dx (k,k) = 



d 4 q 



-9 2 /(V) 



G D (k,k) + ^ 



1 



+ i 

-fc 2 /(V) 



^f x {2TT) 2 5\k + q) 



1 



+ ■ 



+ ^f x (2TT) 2 5\k + q) 



-fc 2 /(V) 



WW} 



eucl. 



(5.10) 



where p = k — k ^ is the total momentum of the qq pair. Returning to Minkowski space 
according to ( |3.5| ) we find the additional contribution to the qq propagator 



G x (k,k) 



i'f; 



u 2 p 2 



l{eW) H [!] 



3 fc 2 /(V) 



mm} 



(5.11) 



and to the axial vector vertex Q4.23] ), 



H^fe, A;) — i 



f: 



{> 



,fe 2 /(V) _ J 2 /(^ 



■>)} . 



(5.12) 

The new contributions 5 X and T^ x to the dressed propagator and vertex are readily seen 
to satisfy the axial Ward identity (4.22). Similarly to the propagator (^^) also the axial 
vertex increases exponentially for k 2 ,k 2 — ► oo, and thus is well behaved only in euclidean 
space. 

As an application of the XSB contribution we evaluate the (AP) correlator between 
the axial vector j§ (x) = q(x)j p J5q(x) and pseudoscalar currents, 

%M v B (p)= [ d 4 xe^^)(0|T{^(x)q(y) 7 5q(y)}|0) 



(5.13) 

which vanishes in ordinary and dressed perturbation theory when m = and f x = 0. 
The correlator is obtained by contracting the axial vector vertex ( |5.12| ) with 75. The loop 
integral is well-defined after a Wick rotation 



mUp) 



fx 



U 2 P 2 



d 4 k 



Tr [e-W^f 



(2tt)4 

and we obtain in Minkowski metric 



mUp) 



757^75 



P 



- fc2 /(2M 2 ) 75| ^ 75 ] eucl. 



p £ + ic 7r 



^/x 2 



(5.14) 
(5.15) 
(5.16) 



This expression satisfies the Ward identity of the correlator, 

p^(p) = -2(0|q(0)q(0)|0) 

given our above result ( [5.7D for the quark condensate. 

According to fl5.15| ) M§(p) is completely determined by the pion contribution (p 2 = 
pole) in our dressed Born approximation. The pion is in fact the only single hadron state 
which couples both to the pseudoscalar and axial vector vertices 5 . As seen from ( |5.16| ), the 

5 The decay constants f nn = for radial excitations 7r„ of the pion, see, e.g., M. 
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behavior of M${p) at large p 2 is given by the quark condensate. This duality between the 
pion and the quark condensate was noted in (Kj . 

It is instructive to evaluate directly the pion contribution to Ai^(p), 



iMUp) 



d 3 q 



(2tt) 3 2|<?| 



d i xe *(*-v) [0 {x o _ y°)(0|i 5 l '(x)|7r(g))(7r( (? )|q(y) 75 q(y)|0) 
+e(y -x )(0\q(y) l5q (y)\7r( q ))(7r(qM(x)\0)] (5.17) 



Using the standard parametrization of the pion coupling to the axial vector current, 

(o\f 5 (x)Hp)} = ///'./•, ( 

and denoting the pion coupling to the pseudoscalar current by 

(O|q(a075q(a0k(p)> = i^e~ ip ' x 

we find 



MUp) = ^—f^ . 

p z + le 

The pion contribution to the (non-amputated) axial vector vertex 
«Tg(x,xi,x 2 ) = (0|T{j5(x)q(xi)q(a;2)} |0) 



(5.18) 

(5.19) 
(5.20) 

(5.21) 



is 



d 3 q 



{i<? U) / d xid x 2 e 



„ik-(x\ —x)^ik-{x—X2) 



(5.22) 



(5.23) 



(2vr) 3 2|q| 

~9{xl - x°)6(x° 2 - x°)e^ x (0\T{ q (x 1 )q(x 2 )} \n(q)) 
-6(x° - x?)0(x° - x° 2 )e~^ x (7r(q)\T {q^^q^)} |0)" 

and involves the non-amputated pion vertex function V w (q,xi — x 2 ) 

(0|T {qafafaixt)} \n(q)) = e'^+^^U W*(q, *1 ~ **) 
In momentum space we find for the singular (pion pole) contributions at (k — k) 2 = 0, 

■Jk^WTT^-^ v ^ 

(5.24) 

V n (k, k) = J d 4 (xx - X2 ) e i(*+*)-(*x-*a)/2 v ^ qf Xl _ X2) (5 _ 25) 
for k — k = q. Comparing ( |5.24 ) with ( 5.12j ) we have 

V v {k,k) 



UKk) = J d A x l d A x 2 ^ k < x ^- i ^- x )Y^{x,x 1 ,x 2 ) 



where 



*Ll. (J- 2 (V) + e *7(V) 



F/ir 



(5.26) 



Thus the 7rqq form factor is exponentially suppressed for k 2 , k 2 



-oo. 
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The pion wave function at the origin ( |5,19| ) is 



^ = -i(0|q(0) 75 q(0)|7r(p)) = -4 J -^ Vn (k,k) 

-^/(W(- IW>+ - P/<V> )-^ 2 - w 

Using this value in ( |5,20| ) we see that the single pion intermediate state approximation of 
the (AP) correlator, M§ w (p), agrees with the dressed perturbative approximation (|5,15| ), 
which was evaluated through a loop integral. 

The pion intermediate state contribution to the (AA) correlator 6 

iDf(p) = J d A xe^ x ^(0\T{j^x)f 5 (y)} |0) = iU 5 (p 2 ) tfgT-W) (5-28) 
is similarly found to be 

n 5 .(p 2 ) = — A- (5.29) 
p z + ic 

and is transverse due to the choice of covariant time ordering. 

The dressed perturbative approximation of the axial vector correlator (AA) ( |5.28|) is 
obtained by contracting the Dirac indices of the qq propagator G(k, k — p) with 7^75 and 
7^75 (respectively) and integrating over the loop momentum k. In Minkowski metric and 
for m = 0, the dressed qq propagator is obtained by adding G x (k, k — p) of (|5.11| ) to ( 4.16j ) 
using Q4.26I) , 

G(k, k) = \[k l - /s 2 p k ] [jg - n 2 p k ] Gj{k\ k\p 2 ) + G x (k, k - p) (5.30) 

poo .4 p —t/i 

"'"•■^-l <5 ' 31) 

Since the Dirac trace of the chiral odd contribution G x (k, k — p) vanishes it does not con- 
tribute to (AA). Using ( [5.28 ) we obtain for Ii§(j> 2 ) the dimensionally regularized expression 

iU 5 (p 2 ) = ^ TT / ^JL(fc-^).(fc-^ )Gfj(fc 3 p (5 . 32) 
J (27r) w 

where A is a renormalization scale. The factor in front of Gi can be expressed in terms of 
the operator Ok = (k — ^ 2 d k ) 2 defined in ( 2.12 ), 

(k - fi 2 d k ) • (k - v 2 d k ) = {k- fi 2 d k ) • (k - fi 2 d k ) = \(O k + 0- k -p 2 ) (5.33) 

Using the identity Q4.19| ) in the expression ( [4.17| ) for Gj we have 

O k G I (k 2 ,k 2 , P 2 ) = ^ (5.34) 



6 The time ordering is understood to be the covariant one described in, e.g., Sect. 5-1-7 of IllJ]. We 
thank the Referee for bringing this point to our attention. 
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This contribution (and the analogous one for O^) vanishes in the dimensionally regulated 
loop integral ( 5,32; ), leaving 



iU 5 (p 2 



3/i 4 



-A 



4-D 



d D k 
(2^p 



(5.35) 



The expression ( |5.31 ) for Gj has p 2 ,k 2 and A; 2 -dependent discontinuities which must 
be specified in the loop integral over k. We shall here only be concerned with the leading 
behavior of IL^(p 2 ) for p 2 — ► — oo, and Wick rotate the integral to euclidean space. Changing 
the integration variable t — > 2uk 2 / u 2 we then find 



d D k 
(2tt) d 



00 exp [-uk 2 /{2fi 2 )] 



(uk - k)< 



eucl. 



(5.36) 



At large (euclidean) p 2 the loop integral is dominated by large k 2 , for which the exponential 
function restricts u < fi 2 /k 2 . Hence uk can be neglected in the denominator, giving 



n 5 b 2 



V 4-D 



3 
1 

6^2 



1 



dPk 



(2vr) D fe 2 A; 2 



eucl. 



(5.37) 



with renormalization dependent constants absorbed in the scale A. 

Although they are suppressed 7 by powers of p 2 at large p 2 , there are non-vanishing 
effects of the dressing on H^(p 2 ), as indicated in ( 5.37| ). On the other hand, the unphysical 
qq propagator Go of fl4.9| ) should not affect gauge invariant quantities at any value of 
p 2 , since it is obtained through a (complex) gauge transformation analogous to (|3,3| ). Its 
implicitly regularized integral over k can indeed be done exactly by choosing the loop 
momentum to be k + q, 



j£k_ 



G D {k~k) 



d 4 k 

d 4 q 



d A q 



(2^ 



(2vr) 2 

-? 2 /(V: 



9 2 /(2M 2 ) 



eucl. 



f d 4 k 


"1" 




h- > 


1 (2tt)4 


A 







eucl. 



f d 4 k 


"1" 




i 


1 (2^) 4 


A 







eucl. 



(5.38) 



yielding the bare (// = 0) result. 



6. Discussion 



We studied some properties of the perturbative expansion of abelian gauge theory around 
a vacuum state not being the empty perturbative vacuum but containing (abelian) gluons 
with vanishing momenta. We defined the expansion through the modification (|1.3| ) of the 
on-shell prescription of the gluon propagator. We showed how the effect of this modification 



7 This is expected, since the dressed propagators approach the bare ones at high virtuality. 



can be summed to all orders in the case of the quark and qq propagators. The generalization 
of our method to higher order (loop) diagrams as well as to other Green functions appears 
straightforward since the modification is equivalent to calculating Green functions in the 
background field fl2,9| ) which is a pure gauge, albeit with an imaginary gauge parameter. 
The complex nature of the field reflects the presence of physical gluons in the vacuum state 
and implies that the dressed Green functions are not equivalent to the standard ones. 

Our modified Green functions may be useful for understanding general properties of a 
confining theory like QCD. They reduce to the standard perturbative result in the short- 
distance limit, but have essentially different (exponential) behavior at long distance. Our 
results on the short-distance limit are already instructive. The safe limit to take is the one 
where all invariants are independently large (and spacelike). The expression ( 4.26| ) for the 



(Dirac reduced, massless) qq propagator indeed turns into a product of free propagators 
when k 2 , k 2 and p 2 are large. However, if we take k 2 — ► oo keeping k 2 and p 2 fixed the 
result is the free quark propagator So(k) multiplied by a long-distance antiquark propagator 
which is distinct from (|2.19| ), which we derived in the absence of a quark. Apparently the 
presence of the quark - even though nearly pointlike - affects the long-distance propagation 
of the antiquark. 

The dressed quark propagator (2.11) has no singularities for p 2 > 0. Thus the pole of 



the free propagator at p 2 = m 2 is removed no matter how small is the parameter A in the 
modification of the gluon propagator ( |1.3| ). This is understandable since the pole arises 
from an infinite propagation distance, during which even a weak background will matter. 
Thus the quark is not an asymptotic (\in) or \out)) state in the presence of a background. 



The discontinuity ( 2.16 ) of its propagator vanishes exponentially for p 2 — > — oo and is thus 
only relevant in the "confinement" region where the quark cannot be treated as a free 
particle. 

Our assumed background field correlates the quark and antiquark giving their joint 
propagator ( f4.27| ) an interesting structure. The kinematic function X(k 2 ,k 2 ,p 2 ) vanishes 
at threshold, y/p* = Vk 2 + vP, where k ■ k = — V k 2 k 2 < 0. Then both terms in 
the denominator of the integrand in ( [4.27 ) vanish within the integration range. As p 2 



approaches threshold from below, this singularity leads to a strong potential between the 
quark and antiquark. Their dynamics becomes effectively non-relativistic since the relative 
3- momentum (see ( |4.28| )) is small compared to the "constituent masses" given by their 
virtualities Vk 2 and Vk 2 . 

It will evidently be important for the usefulness of the present study to identify the 
asymptotic states ("hadrons") and to verify whether the S-matrix formed by them is 
unitary. In section 5 we took a first step in this direction by noting that there is a quark 
propagator solution which breaks chiral symmetry spontaneously in the case of vanishing 
quark mass, and studied some properties of the massless Goldstone boson. 
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A. Connection to scalar abelian gauge theory 



Above we showed that the Dirac structure of the different solutions for the quark propagator 
S of ( 2.11 ) and Sd of ( 3.11 ) could be extracted into separate differential operators that 
operate on scalar functions {Si and Sdi, respectively, see 



jj) below). From Q, (gj) 
we see that in euclidean space Sdi is obtained from the free scalar propagator by Gaussian 
smearing whereas Sd is the smeared quark propagator. This suggests that Si and Sdi are 
solutions for the dressed scalar propagator in scalar abelian gauge theory. In this appendix 
we show that this is the case by checking that Si (and Sdi) satisfy the corresponding 
equation of motion. 

The equation of motion for the scalar propagator Si can be derived similarly to 
Eq. (2.7). The Dyson-Schwinger equation for the dressed scalar propagator is shown in 
Fig. 5. The operator —gdp/fi of Eq. (2J3) may be used to write the equation in an analytic 
form also in the scalar case. As 



1 



/i p p 2 — m 2 fi 2 p 2 — m 2 



- 9 -d° 



-2ig/jLjf 



2ig z g 



2 n a P 



-2igw a ) 



1 



/x 2 p 2 — in 2 



(A.l) 



all gluon insertions are generated: the first line of ( |A.1| ) represents an insertion to a scalar 
propagator and the second line gives an insertion to a one-gluon vertex. The analytic form 
of the Dyson-Schwinger equation is thus 



Si(p) 



1 



/i 2 p 2 

1 1 



+ 



1 



/i 2 p 2 



-2ig\ip a ) (-A 2 



2 n 2 p 2 
1 1 



m- 



■{2ig 2 g a0 ) 



m 



2 (2ig 2 g af s) 



2/i 2 p 2 

The extra factors of 1/2 appearing in 



the last two terms of (A. 2) are either sym- 
metry factors or required to avoid double 
counting of diagrams caused by the second 



order differential operator. (A. 2) simpli- 
fies to 



-A 



2\ 2 



.1? 
ji p 



A- 2 )g a PSi(p) 




- 9 -d^ Si{p) 
9 -d^j Si{p) 



(A-2) 



r~ \ 



[(p _ dp ) 2 - m 2 } Si{p) = O p Si{p) = — 

{A.3) 



Figure 5: The Dyson-Schwinger equation (A. 2) 
for the dressed scalar propagator. 
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Thus using ( |2.13| ) it is easy to check that the scalar part in the quark propag ator (ETUI ), 



i f°° 
s i(p) = 7T~2 / dtexp 
V Jo 



- 2 \v 



m 



l+t 



(A.4) 



indeed satisfies (|A.3| ). The discussion of section 2 shows that it has the correct short 
distance behavior. Similarly, one may check that Gj of ( [4.17 ) is the dressed two-scalar 
propagator. 

It is possible to derive an explicit expression for the perturbation series of Fig. 5. 
Restoring the fi dependence, (A.4) reads 



Slip) 



dtexp 



dtexp 



m 



P 



1 + 



(p 2 — m 2 ) 



exp 



/x 2 i 2 m 2 
2(1 + ^ 2 ) 



(A.5) 



where we scaled t — > t[i 2 . Expanding the integrand at [i 2 = and doing the integration we 

4im 2 fi 2 t 24im 2 (m 2 + p 2 ) fi 4 

(A.6) 



find 

Slip) 



p 2 — m 2 (p 2 — m 2 ) 3 (p 2 — m 2 ) 5 

{-2) k k\ ffJ?\ k j2 



p2 _ m 2 



o oo 

m \ - 

1 + ~^^JT 



2 Ar,2\2fc 



k=l 



m 2 /p 



1=0 



k + l\ (k - 1 
£ + 1 



mr 



p* 



The coefficients in ( |A.6| ) behave as ~ fc!, whence the series is asymptotic (unless m = 0). 
This is expected, as Si has a branch point at /x 2 = (see ( |2.16|) ). 
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